Size and dynamics of vortex dipoles in dilute Bose— Einstein condensates 
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Recently, Freilich et al. [Science 329, 1182 (2010)] experimentally discovered stationary states 
of vortex dipoles, pairs of vortices of opposite circulation, in dilute Bose-Einstein condensates. 
To explain their observations, we perform simulations based on the Gross-Pitaevskii equation and 
obtain excellent quantitative agreement on the size of the stationary dipole. We also investigate how 
their imaging method, in which atoms are repeatedly extracted from a single condensate, affects 
the vortex dynamics. We find that it mainly induces isotropic size oscillations of the condensate 
without otherwise disturbing the vortex trajectories. Thus, the imaging technique appears to be a 
promising tool for studying real-time superfluid dynamics. 
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Introduction — Quantized vortices are topological de- 
fects that carry angular momentum and often serve as 
conclusive evidence for superfluidity. A vortex dipole, on 
the other hand, consists of a pair of vortices of oppo- 
site circulation and can be considered as a basic topo- 
logical structure that carries linear momentum. Al- 
though vortex dipoles are widespread in classical fluid 
flows their role in superfluids seems less well estab- 
lished. Given the appearance of vortices and antivor- 
tices in the Berezinskii-Kosterlitz-Thouless transition Q 
and superfluid turbulence a detailed study of vortex 
dipoles will contribute to a broader understanding of su- 
perfluid phenomena. 

Dilute Bose-Einstein condensates (BECs) are an ideal 
choice for investigating superfluid vortex dynamics 
since they can be described by tractable theories [f| and 
are highly controllable in experiments Q . Vortex dipoles 
have been addressed in various theoretical studies of di- 
lute BECs @-[Hj], and stationary vortex dipoles lol- lil] 
have been found to exist over a wide range of particle 
numbers and trap geometries. Other vortex cluster con- 
figurations that are stationary in nonrotated BECs, such 
as the vortex tripole and vortex quadrupole, have also 
been predicted [Tol - Fl5l ] . In recent experiments 0, [2(J , 
the formation and dynamics of vortex dipoles were in- 
vestigated in 87 Rb condensates. Vortex tripoles have also 
been observed experimentally 21 1. 

Typically, vortices in dilute BECs are imaged by re- 
moving the trapping potential and expanding the con- 
densate in order to make the vortex cores optically re- 
solvable. However, the expansion prevents taking mul- 
tiple images from a single BEC. In Rcf. [2fj| ]. Freilich et 
al. present a method of vortex probing that allows vor- 
tex dynamics to be observed in real time within a single 
BEC rather than by reconstruction of images from multi- 
ple BECs. The method is based on repeated extraction, 
expansion, and imaging of small fractions of the BEC. In 
the demonstration of the technique, vortices were spon- 
taneously created during evaporative cooling due to the 
Kibble-Zurek mechanism, and the orbital dynamics of 



single vortices and vortex dipoles were observed with high 
precision. In particular, clear experimental evidence for 
stationary vortex dipoles was obtained. 

In this Rapid Communication, we present numerical 
results based on the zero-temperature Gross-Pitaevskii 
equation to explain the experimental observations of 
Freilich et al. [20| ■ Previously published predictions 
flit ] for the properties of stationary vortex dipoles concern 
only very weakly interacting or two-dimensional systems 
and hence cannot be used for direct comparison with the 
experiment. Here, we present stationary vortex-dipole 
states for a BEC identical to the one used in Ref. [2(J 
and find excellent quantitative agreement with the ex- 
periment. In order to analyze the accuracy of the multi- 
shot imaging method, we also study the effect of the re- 
peated removal of particles on the dynamics of the vortex 
dipoles, and observe that it excites mainly the breathing 
mode of the BEC but does not significantly disturb the 
vortex trajectories. 

Theory and methods — In the zero-temperature limit, 
the dynamics of a dilute BEC is accurately described by 
the Gross-Pitaevskii (GP) equation Q for the conden- 
sate order parameter ^ , 



iW t ^{v,t) = [H + 3l*(r,t)| 2 ] *(r,t), 



(1) 



where the single-particle Hamiltonian % consists in our 
case of the kinetic energy operator and an axisymmetric 
harmonic trapping potential, 
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V 2 + -mul (x 2 + y 2 ) + -mcu 2 z 2 . (2) 



The atom-atom interaction strength g is related to the 
s-wave scattering length a and the number of atoms N 
by g = 4irh 2 aN/m. Here, m is the atomic mass and 
"J is normalized such that J \ty\ 2 d 3 r = 1. Often, it is 
relevant to consider the dynamics only in the xy plane 
and factor out the z dependence of the order param- 
eter. The general z-factorized form can be written as 
\I/(r, t) = tf)(x,y,t)((z) exp (—ie z t/H), where the function 
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£ is required to satisfy J \(\ 2 dz = 1 and e 2 is a constant 
energy shift. By inserting this into Eq. (Q}, multiplying 
by C* exp (ie z t/h), and integrating over z, we obtain the 
effectively two-dimensional (2D) GP equation 
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where ^2D = 9 J \C\ dz is the reduced interaction strength 
and ^ is normalized according to J \ip\ d r = 1. The en- 
ergy shift is set to e 2 = / (mw 2 z 2 |C| 2 — ft 2 C*C"/ m ) dz/2. 
Although the z dependence of the order parameter does 
not separate ideally in the configuration of Ref. [2(J , this 
technique still provides a useful model for vortex motion, 
as shown in the following. 

Stationary states of the system satisfy the time- 
independent GP equation, which is obtained from Eq. (fTJ) 
or Eq. ([3]) with the replacement ihdt — > fJ-, where /i is 
the chemical potential. For sufficiently large N, a good 
approximation for the lowest-energy stationary state may 
be obtained by neglecting the kinetic-energy term in the 
GP equation. In Eq. (0, this Thomas-Fermi (TF) ap- 
proximation yields the solution 



ii TF (x, y) = Re \ ^n [1 - {x 2 + y 2 ) / R, 



tf] 
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where uq = moj 2 /ngzo is the peak density and the TF 
radius of the 2D BEC is given by Rtf = \J 1 kg-m / itrmJ* . 

In Ref. [2(J , the dynamics of a BEC was observed by re- 
peatedly extracting, expanding, and imaging small frac- 
tions of the trapped condensate across its entire spatial 
extent. The atoms were extracted by a brief (1-3 fis) 
microwave pulse that transferred them to a magneti- 
cally untrapped state, creating a representation of the 
trapped cloud at that instant. We assume that the mi- 
crowave pulse samples the BEC uniformly throughout 
its entire volume, extracting a fraction p of the atoms 
at each spatial point, and neglect possible interactions 
between trapped and untrapped atoms. Under these cir- 
cumstances, the effect of the pulse on the trapped con- 
densate is to reduce the number of particles N by the 
fraction p. In Eq. (j3|), this implies that the interaction 
parameter $2D becomes time dependent, 
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where O is the Heaviside step function and {tj} are the 
time instants of the microwave pulses. 

In the simulations, we search for stationary vortex- 
dipole states in three dimensions and study the dynam- 
ics of symmetric and asymmetric vortex dipoles in the 
effectively 2D limit when particles are gradually removed 
from the trap. The stationary vortex dipoles are solved 
from the three-dimensional timc-indcpcndcnt GP equa- 
tion by using a relaxation method, and the time-evolution 



is solved numerically from Eq. (|3|) by using a com- 
bination of operator splitting and the Crank-Nicolson 
scheme. To obtain dimensionless quantities, we measure 
length in units of the radial harmonic oscillator length 
a r = y/h/mu} r , energy in units of fiw r , and time in units 
of l/u) r . This results in dimensionless forms of Eqs. (U} 
and (|3|) , where the interaction strengths are replaced with 



(3) g = g/hjJrctf. = 4irNa/a r and <?2D = .92D/fe<-VC 
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choose the simulation parameters according to Ref. , 
setting (uj r ,u! z ) = 2ir x (35.8,101.2) Hz and using the 
values a = 5.31 nm [22| and m = 1.44 x 10~ 25 kg cor- 
responding to 87 Rb. Thus, we have a r = 1.80 (im and 
g = 0.0370 x N. 

Results — Figurc[T]summarizes our results for the three- 
dimensional stationary vortex dipolc. For all the parti- 
cle numbers considered, 3 x 10 5 < N < 7 x 10 5 , the 
background condensate density has a parabolic TF pro- 
file along each coordinate direction. The cores of the 
vortex and antivortex are observed to be straight, al- 
though increasing in size toward the surface of the cloud 
[Fig.QJa)]. In Figs.QJc)-QJe), we plot the separation dis- 
tance d sep between the vortex cores, the TF radius -Rtf, 
and their ratio d scp /i?TF as functions of N. Here, d scp 
and i?TF are determined by fitting the z-intcgrated den- 
sity n- m t(x,y) = N J \ty(x, y, z)\ 2 dz to a five-parameter 
surface consisting of a sum of the parabolic TF profile 
[Eq. (gj] and two Gaussians; cf. Fig. [TJb). The separa- 
tion distance d sop is observed to change slowly, decreasing 
by only 2% with N increasing by 133%. On the other 
hand, the TF radius Rtf increases by 16%, thereby be- 
ing mostly responsible for the 18 % reduction in the ratio 
d sep /Rtf- Figure QJe) shows that our numerical result 
for d sep /Rtf is in excellent quantitative agreement with 
the experimental data of Ref. [2(J ■ 

We have also modeled the experiments on the motion 
of symmetric and asymmetric vortex dipoles. For sim- 
plicity, and since no evidence of bent or tilted vortex 
lines was observed in Ref. [see also Fig.QTa)], we have 
computed the dynamics in two dimensions using Eq. ([3]) . 
The value of the reduced interaction strength <72D is cho- 
sen to yield the correct value for the TF radius Rtf- 
The three-dimensional simulation yields Rtf — 15.6 /zm 
at N — 6 x 10 5 , from which we obtain <? 2 d = 4411. We 
have tested the validity of this choice by studying the 
precession of an off-axis vortex located initially at a ra- 
dius r = 0.362i?TF when 5 % of the particles are removed 
every 90 ms starting from t = ms, and found that the 
precession frequency increases from / = 3.15 Hz at t = 
ms to / = 3.61 Hz at t = 655 ms, agreeing well with the 
values / = 3.17 Hz and / = 3.67 Hz reported by Freilich 
et al. [20j. The slight discrepancy in / can be partly 
due to errors in determining the particle number in the 
experiments. 

In Ref. (2ol |. a symmetric vortex dipole was observed 
to remain stationary for at least 655 ms, during which 
nine images were extracted from the BEC. To model this 
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Figure 1: (Color online) Stationary vortex dipole in three dimensions, (a) Isosurface of the particle density |^(r)| 2 for the 
stationary vortex dipole corresponding to the particle number N = 5 x 10 J . (b) Solid blue curve shows the z- integrated density 
along the x axis, 7ii n t(x,0), for the state in (a), whereas dashed red curve is a fit consisting of a sum of the TF profile [Eq. Q] 
and two Gaussians. Here, d sep is the separation distance between the vortex cores and Rtf is the TF radius, (c)-(e) d sep , Rtf, 
and d sop /.RTF as functions of the particle number N for the stationary vortex dipoles. Crosses in (e) represent the experimental 
data of Ref. H. 



experiment, we have computed the evolution of an ini- 
tially stationary vortex dipole when 5 % of the atoms are 
removed from the trap at 90 ms intervals. The results 
are presented in Fig. [21 where Rtf, d sep , and d sep /i?TF 
are plotted as functions of time. In addition to gradually 
reducing the overall size of the BEC, the atom removal 
excites the so-called breathing mode (BM) , causing Rtf 
and d sop to oscillate with the characteristic frequency 

TF 



ip /R 



/bm = 2u! r /2n [Fig. H[a)]. However, the ratio d s 
changes by only 1 % during the first 540 ms. 

Freilich et al. [2(| also studied the precession of asym- 
metric vortex dipoles, with one vortex near the center of 
the BEC and the other close to the surface. We have 
performed a related simulation by first placing the vor- 
tex cores at the measured initial positions and then let- 
ting the state evolve according to Eq. (|3|), with g 2 o be- 
ing reduced by the removal of atoms at regular inter- 
vals [Eq. {5}]. The obtained trajectories are presented 
in Fig. [3] together with the corresponding experimental 
data of Ref. [20[. In Fig. Elb), we also show the tra- 
jectories without the particle removal. We find that the 
vortices move along approximately periodic trajectories. 
The atom removal speeds up the vortex motion and ex- 
cites the breathing mode but does not alter the under- 
lying structure of the paths; it also conserves the orbital 
angular momentum per atom. The outer (anti)vortex 
follows a roughly circular trajectory, whereas the inner 
vortex has a more complicated path with approximate 
discrete rotational symmetry that is sensitive to the ini- 
tial vortex positions; see also Ref. 14j . The separation 




distance d sep between the vortices oscillates strongly with 
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Figure 2: (Color online) Evolution of an initially stationary 
vortex dipole when 5 % of the atoms are extracted from the 
BEC at the time instants indicated on the horizontal axis. 
Panel (a) shows the Thomas-Fermi radius _Rtf (top blue 
curve) and the distance d S e P between the vortices (bottom 
red curve), whereas panel (b) shows their ratio d aep /RTF- 



a period several times shorter than the overall period of 
precession, which in turn causes the instantaneous speeds 
of the vortex cores to vary greatly during the evolution. 

Discussion — In summary, we have studied the size and 
dynamics of vortex dipoles in trapped BECs. Contrary 
to earlier theoretical studies (7-18|, which have consid- 
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Figure 3: (Color online) Trajectories (x v / Rtf , Vv / Rtf) of 
vortices for two different asymmetric vortex dipoles when 5 % 
of the atoms are removed every (a) 90 ms or (b) 60 ms from the 
state (solid blue curves). Dashed red curve in (b) shows the 
trajectories when no particles are removed. Insets depict the 
density \ip(x,y)\ 2 at each moment of particle extraction; these 
time instants are indicated by solid circles on the curves. Tri- 
angles designate the experimental measurements of Ref. [2(| 
for the same time instants and initial configurations. Marker 
colors correspond to the image frames. 



ered only pancake-shaped condensates or small particle 
numbers, we found stationary vortex dipoles in a three- 
dimensional configuration in the TF limit corresponding 
to large atomic clouds. In addition, we investigated the 
evolution of nonstationary vortex dipoles, and found that 
the vortices follow complicated trajectories whose exact 
shapes depend strongly on the initial vortex positions. 

This study explains the recent experiment by Freilich 
et al. [2(j on the real-time dynamics of vortex dipoles. 
Our numerical results for the size of stationary vortex 



dipoles are in excellent agreement with their experimen- 
tal data. We also studied how the atom loss due to the ex- 
perimental imaging method affects the vortex dynamics: 
The major effect is the excitation of the breathing mode 
of the whole BEC, which can be subtracted from the data 
at will by a global distance scaling. Otherwise, the vortex 
trajectories are altered only slightly. Although in some 
cases there may exist more subtle issues such as the onset 
of dynamical instabilities, the multishot imaging method 
used by Freilich et al. appears to be a promising tool for 
studying vortex dynamics in trapped BECs. 
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viding the measurement data and details of the experi- 
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